Demixing in a single-peak distributed polydisperse mixture of hard spheres 
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An analytic derivation of the spinodal of a polydisperse mixture is presented. It holds for fluids 
whose excess free energy can be accurately described by a function of a few moments of the size 
distribution. It is shown that one such mixture of hard spheres in the Percus-Yevick approximation 
never demixes, despite its size distribution. In the Boublik-Mansoori-Carnahan-Starling-Leland ap- 
proximation, though, it demixes for a sufficiently wide log-normal size distribution. The importance 
of this result is twofold: first, this distribution is unimodal, and yet it phase separates; and sec- 
ond, log-normal size distributions appear in many experimental contexts. The same phenomenon is 
shown to occur for the fluid of parallel hard cubes. 
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Our knowledge of the phase behavior of mixtures has 
increased a lot in the last decade. Generally speaking, 
above a certain concentration binary mixtures phase sep- 
arate due to the socalled depletion effect Q . This mech- 
anism has an entropic origin: mixtures phase separate 
when this produces a sufficient gain in free volume to in- 
crease the entropy of the system. More intuitively, deple- 
tion is an effective attraction arising from the unbalance 
pressure resulting when two solute particles are so close 
together that no solvent particle fits in between them. 
Either way we look upon it, in real colloids this mecha- 
nism is supplemented by energetic attractions or repul- 
sions that enhance or inhibit the transition. This is the 
reason why most theoretical research has focused in addi- 
tive hard particle mixtures, in order to ascertain to which 
extent depletion is the basic mechanism of demixing 

Binary hard spheres (HS), the simplest hard-particle 
mixture, has proven very controversial in settling the 
question. More than thirty years ago Lebowitz and Rowl- 
inson Q showed that according to the Percus-Yevick 
(FY) solution of the Ornstein-Zernike equation binary 
HS never demix; a few years later Vrij ^ extended this 
result to p-component HS mixtures in the same approx- 
imation. However, the FY approximation has a well 
known thcrmodynamical inconsistency, as it yields two 
different equations of state, none of which being very ac- 
curate at large densities. It is also well known that a 
linear interpolation of both produces a rather accurate 
equation of state; in the case of mixtures, the latter is 
known as the BMGSL equation of state In spite of 

its higher accuracy, the authors showed that it predicts 
the stability of any binary mixture. 

The BMGSL equation of state is heuristic in its con- 
struction, so this led Biben and Hansen to consider 
the more accurate Rogers- Young closure approximation. 
In this way they found signs of a spinodal instability for 
diameter ratios of the HS smaller than 0.2. Several au- 
thors iQ,^ subsequently confirmed that other approxi- 
mate schemes yield the same instability, although for di- 



ameter ratios strongly dependent on the approximation. 
It was hence believed that phase separation does occur 
in sufhciently asymmetric binary mixtures of HS, an idea 
further supported by the demixing found in binary mix- 
tures of parallel hard cubes (FHG) both in computer sim- 
ulations on a lattice and from fundamental measure 
theory pO[ | in continuum space. 

Presently a different scheme is being accepted to de- 
scribe demixing of very asymmetric binary mixtures of 
hard particles ||ll[]: phase separation occurs between a 
small-particle rich fluid and a large-particle rich crystal 
|l2|-p^ , and this coupling with the translational degrees 
of freedom strongly enhances demixing. It may even 
cause, for very large asymmetry, the appearance of an 
isostructural solid-solid transition as that found 

for narrow and deep attractive potentials — the reason be- 
ing that depletion induces one such potential between the 
large particles ||l5|-|l8t . 

But colloids are hardly mono or bidisperse. By their 
very nature, colloidal particles are usually different from 
each other, and the preparation process of a colloidal sus- 
pension gives rise to a size distribution of particles in the 
colloid. Hence colloids are best modeled as polydisperse 
systems. Polydispersity has received increasing attention 
in the past years ]l9[-p6[| because it causes strong quali- 
tative effects on the phase behavior of the monodisperse 
counterpart. Termination of the freezing transition of 
a HS fluid [^,^, regularization of adhesive potentials 
p5| , p^ , or appearance of a vapor- liquid transition in ad- 
hesive HS are but a few examples. 

Phase equilibrimn in polydisperse systems is more 
complicated than in monodisperse systems: equilibrium 
equations between coexisting phases become functional 
equations when polydispersity is present Much 
of the effort has indeed concentrated in taking advan- 
tage of extra symmetries (as the dependence of the excess 
free energy on the moments of the distribution p3| , p^ ), 
or in mapping the polydisperse system into a simpler one 
(e.g. a binary mixture [p2[), with the purpose of reducing 
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equilibrium to a few algebraic equations. Specifically, by 
means of a moment-based formalism |2^,Q Warren |2^] 
has shown that while a binary fluid of HS in the BMCSL 
approximation never demixes, it does if enough polydis- 
persity is introduced and the diameter ratio of the two 
"main" species is sufficiently small. This unexpected re- 
sult proves that the effect of polydispersity can be more 
subtle and nontrivial than one can tell a priori. 
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FIG. 1. Distribution, p{s) (in percentage per /im), of diam- 
eters, s (in /im), for a suspension of mullite (3AI2O3 ■ 2Si02), 
a ceramic powder used in structural ceramic applications. 
Bullets are the experimental data; the line is a fit of a 
log-normal distribution of mean x = 1.57 ^.va and variance 
w — 0.43^. (Experimental data courtesy of R. Moreno, Insti- 
tuto de Ceramica y Vidrio, CSIC, Madrid, Spain.) 

In this letter I will push this result much further. 
Though polydisperse, Warren's system has two clearly 
differentiated species. I will show that for demixing to 
occur this is not necessary if the size distribution de- 
cays sufficiently slowly. Again using a moment-based 
formalism (equivalent to Warren's but based on different 
grounds) to obtain a spinodal of an arbitrary mixture 
(with several species, polydisperse, or both), I will prove 
that BCSML polydisperse HS with a log-normal size dis- 
tribution (i.e. the logarithm of the diameters follows a 
normal distribution) phase separate provided the vari- 
ance of the distribution is above a threshold. This result 
is surprising because one such distribution is unimodal. 
It is thus clear that the long-size tail plays a crucial role 
in inducing demixing, but although slowly decaying, this 
tail decays faster than any power-law. In addition, it is 
important to stress that, contrary to what one might be- 
lieve, log-normal distributions are very common in many 
polydisperse systems: ultrafine metal particles, paint and 
rubber pigments, photographic emulsions, dust particles, 
aerosols, cloud droplets, or suspensions of ceramic pow- 
ders (see Q and references therein; see also Fig. |l|), etc. 



Therefore this effect should be experimentally accessible, 
provided any of these systems can be prepared with suf- 
ficiently large size dispersion. 

In what follows I will briefly describe the formalism to 
obtain the spinodal in the case that the excess free en- 
ergy of a mixture depends on a few moments of the size 
distribution. Then I will apply this formalism to poly- 
disperse HS and PHC in scaled particle approximation 
(equivalent to PY in the case of HS), and I will finally 
prove the above mentioned result. 

Suppose we have a multicomponent hard particle mix- 
ture, whose free energy per unit volume is / = [3F/V = 
jid ^ jex^ ^i^j^ pA ^ pJlogKp, - 1], V, the thermal 
volume of species i, (3 the reciprocal temperature (in units 
of the Boltzmann constant), and pi the number density 
of species i. Let us further assume that f'^'^ — (f>{£,), 
^ denoting the set of moments {■Coj Cij • ■ • j Cp}; where 

— J2i ^iPi i'^i is the diameter of species i). Stability 
of the mixture requires / to be a convex function in all of 
its variables. As in any hard particle system / depends 
on temperature through an additive term (the free en- 
ergy of the ideal gas) and hence it is always convex with 
respect to this variable. The study of stability is thus 
restricted to the set of densities of the species involved. 
If we define the matrix 

M = Tr^--% + — ^, (1) 
opidp-j pi opiopj 

then stability requires it to be positive definite. As it is 
so for low densities (ideal mixtures are always stable) the 
spinodal is usually described as the set of points in the 
densities space where |M| = (| • | denotes the determi- 
nant) This condition is equivalent to M having a 
zero eigenvalue, i.e. M • u = for some vector u 7^ 0. Let 
us define for convenience a new vector e as Mi = piCi; 
then M • u = means 

-pjCj. (2) 
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For a polydisperse mixture, if p{s) = pp{s), where p is 
the total number density of particles and p{s) is the size 
probability distribution (in terms of a dimensionless di- 
ameter s, for instance), /"^ = / ds p{s)[logV (s) p{s) — 1] 
p9[ , ~ J dss''p{s), but f^^ is still assumed to be of 
the form f''^ — </>{£,). The stability condition is now the 
positive definiteness of the integral operator whose kernel 
is defined by 
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(3) 



While the determinant condition has no direct equiva- 
lent for an integral operator, the zero-eigenvalue one is a 
straightforward extension of (|^) : 

sy 



dt 



6p{s)Sp{t) 



p{t)e{t) 



(4) 
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But we can simplify this equation by using 



which transforms Eq. (Q) into 

p 

e{s) = - <^kis''ai, $fc; = 

k,l=0 



(5) 
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where ai = J ds p{s)e{s) . Substitution of Eq. into 
the latter definition yields pll] 



fc,;=o 



(7) 



For this equation to have nonzero solutions for a we must 
have |Q| =0, with 
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mn — ^mn 



'Y^m+k^kn, m,n = 0,...,p. 



k=0 



(8) 



For convenience we introduce the variables yk = 
~ ^3) and define two new matrices, Y and fl, by 
Ym„ = j/m+ri and = (1 - ^3)$m„. Then the condition 
reads T) = \I + Y ■ n\ =0, with I the identity matrix. 

Many scaled-particle theories can be written 0,|^,|2| 
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(Notice that if (j) and are all multiplied by the same 
constant the product Y • fl remains invariant, so we will 
assume these variables defined up to a constant.) For 
this particular choice of we will have 



n = 



fo 1 

A Ay2 

A 6By2 Ayi+&Byl 

\ 1 Ay2 Ay I + QByl yo + 2Ayiy2 + &Byl 
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It is then straightforward to check that 



[1 + (1 - A)^,f + {&B A^)i2U 

(1-6)^ ' ^ ^ 

where we have replaced the y's in terms of the ^s; ^ is 
the volume fraction of the fluid. We are now ready to 
discuss specific cases. 

Percus-Yevick HS correspond to choosing A — ?>, 
B — i/2 (definitions of (j) and ^ carry an extra 7r/6 fac- 
tor); then V = (1 + 2^)2/(1-^)4^ the resuh obtained by 
Lebowitz and Rowlinson |2| for the binary mixture and 
later generalised by Vrij [I for a multicomponent mix- 
ture. The validity of this result has now been extended 
for a polydisperse mixture [Bll. 



Scaled-particle free energy for PHC is given by [|lO 30 
A = 3, B = 1. Then 



2? = [(1 + 2e3)' - 36^41/(1 - 6)' 



(12) 



By a suitable choice of ^ and ^ we can make 2? = 0. If 
we consider a binary mixture we recover the result of Ref. 
|l0| . By defining mk = £,k/(,3— J dss^p{s), the mixture 
of PHC in this approximation will be stable provided 
2^(6) > m2m4, where 7^(6) = (1 + 26)7(3^1)- The 
function K{x) diverges at x = and it monotonically 
decreases to 3 as ^ approaches 1. Then 17121714 > 3 is the 
condition to find demixing, and the larger the product 
7712^4 the smaller the value of ^ at which it appears. 

Schulz distribution is a common choice in studying 
polydispersity ||l^. If the mean is set to 1, it can be 
written as 



p{s) = w-^T{w-^){w''^sy 



(13) 



with < w < 1 the variance. Its moments are ^„ — 

pYlk=o(^ + ^''^'^) fo'^ - 1' m2m4 = (1 + 3w2)/(l + 
2w'^) < 4/3. Therefore there is no demixing for such a 
distribution. It is easy to check that there is no demixing 
either for a uniform distribution in any interval of diam- 
eters. However if we consider a log-normal distribution 
with mean 1 and variance < w < cxd. 



V27rln(l + ^2) 



exp 



ln'[(l 



21n(l + w2) 



(14) 



whose moments are given by = ^(l-l-w^)"^"^^)/^^ then 
7712^4 = 1 + w'^ , which, by increasing w, can be made 
arbitrarily large. There is thus demixing for a log-normal 
distribution provided w > \/2. 

The latter result rises the question whether the same 
occurs for BMCSL HS. In this case the free energy is 
given by pi 



72 ^ ^0 
?3 



ln(l - 6) + 



?2 



366 
1-6 6(1-6) 



(15) 



where, as in PY, (j) and 6 are defined with an extra 7r/6 
factor. The difference with respect to the PY one is 



= el 



^ln(l 



6) 



1 - 36/2 

6(1-6)2 



(16) 



where = 0bmcsl — 4>'py- The corresponding matrix 
f^BMCSL — f^PY + Ail, where Afi is obtained from the 
second derivatives of A0 (its only nonzero elements Ailj-fe 
are those with j, /c = 2, 3). 

The resulting expression for T) with the moments re- 
placed by those of a log-normal distribution is a com- 
plicated formula relating the volume fraction, 6? and 
the variance, w\ nevertheless, the spinodal can be nu- 
merically determined (Fig. 0). It can be seen that the 
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threshold variance for demixing is w « 1.6. The inset 
shows the size distribution for a variance for which the 
spinodal occurs at hquid densities, in order to illustrate 
the long tail at large particle sizes. Comparing with Fig. 
|l| and with the figures of Ref. Q we can see that stan- 
dard experimental samples are still far from reaching the 
spinodal. 
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FIG. 2. Volume fraction, ^3, as a function of the variance 
of the log-normal size distribution, w, at which a spinodal 
instability appears. Inset: size distribution for the particular 
value w — 2.5. 

In summary, I have introduced a simple formalism to 
compute the spinodal of fluids whose excess free energy 
can be described as a function of a few moments of the 
size distribution of a mixture. By means of this formal- 
ism I have proven that a polydisperse mixture of PY HS 
never demixes despite its size distribution. I have also 
shown that a polydisperse mixture of PHC may demix 
only if the size distribution is long-tailed at large sizes, 
as for instance a log-normal distribution. More impor- 
tantly, polydisperse BMCSL HS, also log-normally dis- 
tributed, do demix for a sufficiently large variance. Re- 
markably, log-normal distributions have a single charac- 
teristic size, so the driving mechanism of the transition 
in these systems is not as clearcut as in the case of binary 
mixtures. Finally, sizes in many polydisperse systems are 
log-normally distributed. Present experimentally avail- 
able systems, though, are far from the spinodal as pre- 
dicted by this work. I hope that this result encourages 
experimental work to achieve this limit and verify the 
conclusions reported here. 
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